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CALCULATION OF WIND-DRIVEN SURFACE CURRENTS 
IN THE NORTH ATLANTIC OCEAN 

Thomas H. Rees and Richard E. Turner 
Langley Research Center 

SUMMARY 

Calculations to simulate the wind-driven near-surface currents of the North Atlantic 
Ocean during Northern Hemisphere summer are described. A computer simulation pro- 
gram, based upon a primitive equation formulation, was used to integrate the equations of 
motion on a finite -difference grid with 2.5° resolution. The model ocean was homogeneous 
with a uniform depth of 100 m, with five levels in the vertical direction, and with a con- 
strained upper surface. The equations were integrated for approximately 40 simulation 
days, at which time the computed surface currents agreed generally with observed cur- 
rents. The development of a subsurface equatorial countercurrent at the 90-m level was 
observed. 


INTRODUCTION 

A number of problem areas exist today which involve the interaction between the 
environment and man’s activities. Of particular concern are the injection and subsequent 
dispersal of manmade pollutants in the atmosphere and in the surface waters of the globe. 

Consequently, an effort was begun at the Langley Research Center to develop a gen- 
eralized circulation model suitable for simulation of both atmospheric and oceanic circu- 
lation on a wide range of scales. A generalized tensor formulation of the primitive equa- 
tions (conservation of mass and momentum and transport of conservable scalars) was 
derived in reference 1. This formulation provides a set of transport equations common 
to the atmosphere and ocean by using similar coordinate systems and by assuming both 
to be in hydrostatic equilibrium in the vertical direction. Separate but similar computer 
models based upon the generalized formulation were written for the atmosphere and the 
ocean. A detailed discussion of the ocean dynamics model is presented in reference 2. 

Because large-scale changes in the open ocean occur slowly, a common practice in 
ocean modeling to improve computational efficiency is to adopt a "rigid-lid" approximation 
in which the vertical velocity is set identically to zero at the sea surface. This effectively 
removes surface gravity waves, without appreciably affecting the steady circulation pat- 
terns, and thereby allows time increments to be orders of magnitude larger than would 
otherwise have been possible with simple explicit numerical schemes. However, the 



application of a rigid lid has in the past generally led to formulation of the equations of 
motion in terms of a stream function. (See ref. 3, for example.) A consequence of the 
stream function approach is that a Poisson boundary value problem must be solved for the 
stream function around land-mass boundaries. Islands greatly complicate the solution. 

The model of reference 2 combines a general formulation of the primitive equations 
with a new rigid-lid approximation. Land boundaries are handled much more simply than 
in a stream function model, while the advantage of a longer time step is retained. 

The purpose of the present work is to provide preliminary validation of the ocean 
circulation model of reference 2. For the purpose of this report, preliminary validation 
consists of (1) demonstrating long-term stability of the computer model, (2) demonstrating 
convergence of the model to an equilibrium state, and (3) verifying that the results of the 
calculations agree with the observed data to an accuracy commensurate with the model 
spatial resolution. 

The approach taken was to define an idealized model of a closed basin. The cir- 
culation of the North Atlantic Ocean was chosen for simulation because of the availability 
of observational data and because at least the surface layer could be isolated from other 
bodies of water without severely affecting the results. A grid spacing of 2.5° in longitude 
and latitude was chosen as a compromise between finer resolution and increased computer 
requirements. The ocean was idealized to be homogeneous in density, to be shallow and 
of uniform depth, and to be driven only by a summer mean wind stress distribution. 

SYMBOLS 


A 

Cd 

Fa 




K 


H 


Kr 


K 


V 


determinant of quasi-horizontal metric tensor 
drag coefficient for bottom friction calculations 

mixing, centrifugal, and incremental Coriolis force terms, defined by equa- 
tions (4) and (5), N 

Coriolis parameter 

covariant element of the metric tensor 
horizontal subgrid mixing coefficient, m2/sec 
Von Karman constant, 0.4 
vertical subgrid mixing coefficient, m2/sec 
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mixed tensor components of deviatoric strain-rate tensor, defined by 
equation (8) 

pressure, N/m2 

surface constraint pressure, N/m^ 

approximate hydrostatic pressure referenced to mean sea level (flat top) 
geoid, N/m^ 

mean radius of the Earth, m 
salinity, parts per thousand 

physical velocity of coordinate grid points relative to rotating Earth, m/sec 
temperature, °C 
time, sec 

physical component of fluid velocity relative to rotating Earth, m/sec 
turbulent part of Uj, m/sec 
approximate value of Uq,, m/sec 

physical component of fluid velocity relative to coordinate frame of reference, 
m/sec 

depth of ocean in transformed coordinates 
coordinates of reference coordinate system (see fig. 1) 
mixing length for subgrid diffusion 
generalized density, P\js^ 

specifies layer thickness distribution as a function of x3 
density of sea water, kg/m 3 



reference density, 10^ kg/m3 
components of wind stress, N/m^ 
i// variable of integration 

CO angular velocity of the Earth, rad/sec 

Indices: 

i,j take on values 1, 2, or 3 

Of takes on values 1 or 2 only 

1,2,3 particular coordinate direction 

Notation: 

< > space average over a grid cell and time average over a time step 

( ) i partial derivative of ( ) with respect to xi; for example, P o = 

’ 9x2 

MATHEMATICAL MODEL 

The generalized circulation model described in reference 1 is a tensor formulation 
of the so-called primitive equations, the equations representing conservation of mass, 
momentum, salinity, and internal energy. The formulation makes use of the hydrostatic 
approximation, which neglects local acceleration in the equation of motion for vertical 
velocity; the momentum equation in the vertical direction is replaced by the hydrostatic 
and continuity equations. This approximation is quite common in large-scale circulation 
modeling because it allows a much longer integration time step. Some notable features 
of the mathematical formulation are 

(1) Combining density and the vertical scale factor to form a generalized density 

(2) Distortion of the quasi-horizontal bottom coordinate surface to follow the Earth's 

topography and provision for vertical motion of the other quasi-horizontal 
surfaces 

(3) Projection of gravity onto the quasi-horizontal surfaces 

(4) Derivation of an improved subgrid diffusion equation 
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(5) Use of generalized density rather than pressure or density as a dependent 

variable 

The ocean circulation model described in detail in reference 2 is an extension of the 
generalized formulation of reference 1. The following two approximations, both of which 
are common in ocean modeling, are made: 

(1) The Boussinesq approximation, which neglects density variations except in the 

buoyancy term and in the hydrostatic equation 

(2) An empirical equation of state, giving density as a function of pressure, temper- 

ature, and salinity, is assumed. 

The ocean model has the capability to operate as an unfiltered primitive equation model 
(as a free-surface model) or as a filtered model (with the rigid-lid approximation). Land 
boundaries and bottom topography are modeled simply. 

hi this report the dynamic equations are presented in somewhat simplified form. For 
example, since the present study makes use of the rigid-lid approximation, terms dealing 
with motion of the grid structure are omitted. Terms in the momentum equations arising 
from variation in depth have also been omitted. Since the present study deals with a homo- 
geneous ocean, the transport equations and boundary conditions for salinity and internal 
energy are omitted. Equations have been combined when feasible. The reader interested 
in a rigorous derivation of the generalized model is referred to .reference 1, and the reader 
interested in a more detailed description of the ocean model, including the rigid-lid approx- 
imation and the numerical technique employed, is referred to reference 2. 

Coordinate System 

The coordinate system (fig. 1) chosen for the ocean model is described in detail in 
reference 1. The lower coordinate surface (x3 = o) follows the ocean floor, while the 
upper coordinate surface (x3 = X^) follows the upper surface of the ocean. For the pres- 
ent case of a uniform-depth ocean with a constrained upper surface, the reference frame 
reduces to spherical polar coordinates. The coordinate xl represents colatitude and 
x2 represents longitude. The gravity vector coincides with the local x3 parametric 
line. 

In appendix B of reference 1, it is shown that the differential geometry of the coor- 
dinate system is sufficiently represented by the approximate metric tensor 

gll = r2 

g 22 = (r sin xl) 

S33 = 
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where r is the Earth's radius, assumed constant everywhere. The determinant of the 
quasi -horizontal metric tensor is simply 


A - giig22 


The conservation equations are formulated in terms of a generalized density 
and appl3dng the Boussinesq approximation /p = p J results in 




^ = PQ\J^ 


where ^ggg is the vertical scale factor. The rigid-lid approximation described in ref- 
erence 2 negates long-term depth variations, so that the time-averaged motion of the coor- 
dinate system relative to the rotating Earth is zero. However, although the horizontal 
components of grid velocity Sq, are negligible, the instantaneous vertical grid veloc- 
ity S3 must be considered, so that 




u 


a 


V3 = Uj - S3 

The spatial derivative S3 is ignored. 

Dynamic Equations 

If substitutions are made according to the foregoing approximations, then the con- 
tinuity equation given in reference 1 as equation (1) may be written as 


iT 


r sin x-* 




^ Po^3,3 = 0 


( 1 ) 


The second and third terms on the left side of equation ( 1 ) represent horizontal and verti- 
cal divergence, respectively. Similarly, the d3mamic equation for momentum conservation 
in the x^-direction (south direction) given in reference 1 as equation (2) may be written as 




= --L. 


^ P,1 ^ ?fU2 - 


( 2 ) 
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where the Coriolis parameter f is given by 
f = 2(t) cos 

and where the subgrid diffusion and centrifugal force terms are grouped into Fj. The 
corresponding equation in the x2 direction (east direction) given as equation (3) of 
reference 1 may be written as 


9t\ 




p^r sin X 


: P 2 - ?fUi - F2 

1 r» V-*- 


(3) 


The terms Fj^ and F 2 from equations (2) and (3) may be written explicitly as 
1 


Fi = 


r sin x-* 


sinxl)^^-. (?<U’iU^>)^J 

- (U 2 U 2 + <u^u^))L£^ 


(4) 


^2 = 


r sin X-* 


? cot x^ 


+ p. 


,(<«2U3>) 


,3 


+ (U1U2 + <UiU 2 )) 


(5) 


The terms on the right side of equations (4) and (5) represent horizontal and vertical sub- 
grid mixing and the contribution of U 2 to centrifugal and Coriolis forces. The primed 
quantities indicate the turbulent parts (i.e., the difference between the instantaneous values 
and the space-time averages), and the brackets < > indicate space averaging over a nodal 
cell and time averaging over a time step. The horizontal subgrid mixing is parametrized 
by a nonlinear kinematic eddy viscosity theory, wherein the subgrid velocity correlations 
are, from equations (8) and (9) of reference 1, 


<uiul> = 

<“l“2> = 


-KhJ^i 

%®*2 


Sin X-*- 


(6) 
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The kinematic eddy viscosity Kjj given in equation (11) of reference 1 may be expanded 
and rewritten as 




1/2 


(7) 


where Kq is the Von K^man constant for turbulent mixing, and A is a characteristic 
mixing length. The required components of the strain-rate tensor are, from equation (7) 
of reference 1, 


Pi = ^ 


r sin x-* 




sin x^^ ~ ^2 2 ~ ^^1 


f ^ = — 
2 r 


|u 2 sin x^j ‘*'^12“ ^^2 


« 1. 


(sin x^j 
1 


^ 2 = 1 




( 8 ) 


The vertical subgrid mixing is parametrized by a much simpler technique, wherein a 
mixing coefficient is specified. By combination of equations (12), (13), and (7) of 
reference 1, the velocity correlations in the vertical mixing terms of equations (4) 
and (5) may be written as 


(m) = -^“ 2,3 


(9) 


where the mixing coefficient Ky may vary with depth. At the lower boundary, the 
conditions, from reference 2, 



O t 

>ll/2 ^ 

<u;u^ = -CdUi 




O c% 

■1/2 [ 

<U2«3> = -^ d '^2 






are applied to equations (9), while at the upper surface the wind stress 


( 10 ) 


= Po<"l''3) 


^2 = '>o<''2«3> 


is a driving boundary condition (forcing function). 
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( 11 ) 



The vertical velocity Vg required in equations (1), (2), and (3) at a given location 
xl,x2] may be computed as 



= -S3(x 3) + Ug(x3) (12) 

where ^(x^) is the ratio of the local height above the sea floor to the total ocean depth. 
Equation (12) follows from equation (15) of reference 1 with the boundary conditions, 

U 3 =Vg =0 at the lower boundary (x3 = o) 

Vg = 0 at the upper boundary (x3 = x^) 

and with the restriction of equally spaced layers, or more generally, with S3(x3) pro- 
portional to the local height above the sea floor. 

Land-Sea Boundaries 

In the present work, the boundary condition applied at model land-sea boundaries is 
simply that the mass transport normal to such a boundary in a water column adjacent to 
the boundary is zero. This condition is automatically satisfied by requiring the depth to 
slope to zero at a nodal point instead of between nodal points, so that ^/ggg and ? vanish 
identically along the beachline. Thus, the coastline always follows a parametric line con- 
necting two nodal points. The continental boundaries in this system are series of steps. 

When the rigid-lid approximation is used, an additional boundary condition — that the 
mass transport normal to the beach in a water column adjacent to a land boundary is 
zero - is required. This condition is imposed by adjusting the horizontal velocity com- 
ponents throughout the column to yield a net transport of zero for the column. 

Air-Sea Interface 

For the present study, the wind stress on the upper surface was chosen as the driving 
condition. Other boundary conditions at the air-sea interface (evaporation, precipitation, 
and heat transport) were zero. 
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Surface Constraint Pressure 
In reference 2, it is stated that the pressure is given by 


where is an approximate hydrostatic pressure referenced to the mean sea level geoid 
(i.e., the flat top depth) and P^ is the surface constraint pressure. In the empirical 
equation of state, p = p(P,S,T), the compressibility effects of P^, are ignored. Thus, 
for the case of a homogeneous, uniform-depth basin as presented herein, Pf q, = 
in equations (2) and (3), P qi = Pc,o!‘ surface pressure P^ is computed iteratively 
as described in reference 2. 

NUMERICAL SOLUTION 

The mathematical model was programed for a digital computer in finite-difference 
notation and integrated on a three-dimensional grid. This integration technique used 
central differencing in time and space, with the horizontal grid structure split into two 
time layers. 

Application of the rigid-lid constraint is achieved by dividing the integration into 
two distinct stages. The first stage computes approximate values of Ua over the entire 
lattice by setting P^. = 0. In the second stage, these approximate velocities are corrected 
to enforce zero local horizontal divergence. The sequence of computations is summarized 
as follows: 

(1) V3 is computed from equation (12). 

(2) is computed from equation (1). 
at 

(3) are computed from equations (2) to (11) with P^ = 0 (no pressure 
gradient force). 

(4) ? is updated. 

(5) Uq; are updated to approximate values Uq,. 

(6) Boundary conditions at lateral land-sea boundaries are applied. 

(7) Pg is computed by relaxation to yield values of -^(^Uq!) which force the 

solution to zero local horizontal divergence. 

(8) U(x are updated to the corrected values U^. 

The numerical technique is discussed in more detail in reference 2. 


10 



NORTH ATLANTIC BASIN MODEL 


A lattice spacing that would give reasonably good definition of large-scale flow 
features was sought. A 2.5° by 2.5° grid spacing in longitude and latitude was adopted as 
a compromise between finer resolution and increased computer requirements. The 
decision to model the North Atlantic Ocean in isolation from other oceans required the 
use of either open boundaries or artificially closed boundaries at the north and south grid 
limits. The latter were chosen for simplicity. Also, a closed system was considered to 
have advantages for validation purposes. 

The choice of artificial boundary locations was not arbitrary. Boundaries were 
sought at which the flow along the boundary would, in the real world, be predominately 
parallel to the artificial boundary and/or at which the mass transport across the boundary 
would be minimal. It was recognized that near 7° to 8° south latitude the surface current 
(the South Equatorial Current) is essentially zonal. (See fig. 2 which was reproduced from 
ref. 4 by permission of publisher.) This location was therefore deemed a reasonable 
choice for an artificial boundary in a surface-layer model. (Since the deep currents near 
the Equator are predominately meridional, the 7° to 8° latitudes would not have been a 
good choice for a closed boundary for a deep-ocean model.) For the northern boundary, 

80° north latitude was considered a reasonable choice because the north-south surface 
transport is small there. The grid used for the present calculations is shown in figure 3. 
As indicated, the grid encloses the region from 82.5° north to 10° south and from 100° west 
to 20® east. The computational grid was composed of 48 points in the east-west direction 
and 37 points in the north-south direction. 

MODEL PARAMETERS 

The input parameters for the present study are summarized in table I. The lattice 
spacing of 2.5° in longitude and latitude is equivalent near the Equator to a square mesh 
with 278 km between grid points, whereas at 80° north the resolution in the east-west 
direction is approximately 48 km. (This resolution tends to smear the boundaries of 
strong currents such as the Gulf Stream.) The model had five layers for the present study 
with a uniform depth of 100 m. Each layer was 20 m thick. 

For the present calculations, the wind stress field was interpolated from tables of 
the summer mean presented by Hellerman in reference 5. The applied wind stress vector 
field is shown in figure 4. The pattern of figure 4 represents the Northern Hemisphere 
summer mean wind stress based on wind-rose data given in U.S. Navy atlases. The wind 
stress data in reference 5 were tabulated on a 5° grid and extended only to 72.5° north. 

The data presented in figure 4 were interpolated onto the present grid system. Note, 
however, that the data were not extrapolated into the most northern latitudes, and that the 
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applied wind stress east and west of Greenland is very low. (These areas are largely 
covered with ice in the real world.) 

Because of the high level of turbulence in the wind-driven surface layers, the verti- 
cal mixing coefficient for momentum transport Ky was taken as 0.01 m2/sec. This 
value of Ky agrees with that used for the upper 100 m in reference 3. The Von Karman 
constant Kq in equation (7) was taken as 0.4. The drag coefficient Cq in equation (10), 
normally taken as approximately 0.002, was taken herein as 10"^ in an attempt to simulate 
the lower resistance of the water at 100 m as opposed to the resistance of a real bottom. 

The salinity and temperature fields were set to standard ocean conditions = 35 
and T = 0° cj and held constant. The velocity field was initially zero. The kinetic 
energy in the entire field was monitored as a convergence check. 

RESULTS 

The dynamic equations of motion were integrated numerically on a digital computer 
starting from a state of rest. The initial response to the wind stress was development of 
a large, slow-moving clockwise gyre centered near the geographic center of the ocean. As 
simulation time progressed, this gyre strengthened, the gyre center moved westward, and 
the western boundary currents intensified. 

Three horizontal velocity vector fields are presented in figure 5. This figure repre- 
sents the computed, nearly steady-state currents after approximately 40 simulation days. 

The currents in the surface layer (0 to 20 m) are shown in figure 5(a); those in the bottom 
layer (80 to 100 m), in figure 5(b); and the vertically averaged horizontal currents, in 
figure 5(c). The time derivatives from the momentum equations were smoothed zonally to 
control numerical problems arising from convergence of the x^ parametric lines near 
the pole. This smoothing was applied heavily at the most northern latitudes, but was not 
applied at all near the Equator. 

Toward the end of the simulation, significant numerical problems arose in the region 
east of Greenland from about 60° north to the northern limit of the model. Within 3 simu- 
lation days after the time period shown in figure 5, the velocities along the two most north- 
ern nodal lines began increasing exponentially. For approximately 12 days prior to this, » 

however, the entire field was essentially in steady state. The high-latitude instability prob- 
ably occurred because of the minimal applied wind stress in that region (fig. 4). The dis- , 

organized motion in the region was driven by diffusion, and thus no predominant pattern 
existed. 

Two other aspects of the model, grid resolution and land boundary treatment, warrant 
attention. A second-order accurate finite -difference technique such as the one used herein 
tends to spread strong gradients over about five nodal points. Thus, strong currents (the 
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Gulf Stream, for example) are much wider in the computed current field than in the real 
world and have a lower maximum current speed. The land boundary treatment in the 
present model is to set the vertically integrated horizontal transport momentum normal 
to a boundary equal to zero by adjusting the horizontal velocity in each layer. Tests for 
land boundaries are carried out along x^ and x^ parametiic lines, and the model sees 
irregular continental boundaries as series of steps. Consequently, strong boundary cur- 
rents are affected in three ways. First, since the momentum adjacent to the boundary is 
reduced at a boundary step, the momentum of the current may be significantly reduced. 
Second, reducing the current velocity next to a land boundary increases the shoreward cur- 
rent gradient. As a result, on a finite-difference grid, the current is forced out from the 
boundary by approximately one nodal distance. And third, a step change in the lateral 
boundary, especially a double step, can cause a strong boundary current to be deflected 90®. 

The prominent features of computed surface current patterns shown in figure 5(a) 
agree for the most part with the corresponding features of the observed data. (See fig. 2.) 

A rough quantitative comparison of the computed and observed current speeds for some of 
the major currents is presented in table n. 

In figure 5(a), the westerly South Equatorial Current (1) in the most southern lati- 
tudes is well developed, even though in the model it must be fed principally by overturning. 
(The Benguela Current, which flows northward off the coast of South Africa to feed the South 
Equatorial Current in th§ observed data, is not modeled because the southern boundary is 
closed.) The South Equatorial Current splits near the east coast of South America, with 
part turning back to form the Equatorial Countercurrent (2) and part continuing up the 
coast. The weakness of the coastal current may have been due in part to the stepwise 
representation of the land geometry and land boundary treatment which combine to turn the 
flow northward. The observed South Equatorial Current (fig. 2) also splits, the southward 
branch forming the Brazil Current. The wind stress (fig. 4) along the three most southern 
latitudes of the model is directed toward the west- northwest. Since these latitudes (1.25® 
to 6.25® south) are in the Southern Hemisphere, Coriolis acceleration acts to the left; the 
result along the closed boundary is, a current with a small southerly component. The flow 
is blocked by the artificial southern boundary, however, and merely sinks. The North 
Equatorial Current (3) splits also, with part flowing south of Cuba and the rest to the north 
of Cuba. The two segments rejoin in the Florida Current. The Gulf Stream (4) and the 
North Atlantic Current (5) are also easily recognizable. The Subtropical Convergence (6) 
(the gyre center) is well situated and. well defined. 

The computed currents at the 80- to 100-m level, shown in figure 5(b), differ little 
from the surface currents, except in the equatorial region. At the lower level, the South 
Equatorial Current (1) is considerably weaker, whereas the Equatorial Countercurrent (7) 
is more pronounced. The development of a stronger subsurface equatorial countercurrent 
is interesting since such a phenomenon was not expected in a homogeneous model with 
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fairly coarse grid spacing and with a shallow, flat bottom. Furthermore, the model 
parameters tend to oppose development of subsurface countercurrents; the high vertical 
mixing coefficient (o.Ol m2/sec) and the low bottom drag coefficient (lO"^) tend to reduce 
vertical gradients in the velocity field, as can be seen in the rest of the region. 

The vertically averaged currents in figure 5(c) show more clearly the computed cir- 
culation. The dominant features of the flow are three fairly well-defined gyres. The 
large anticyclonic gyre including the Gulf Stream, the North Equatorial Current, and the 
North Atlantic Current is bounded on the south by a smaller anticyclonic gyre, and on the 
north by a cyclonic gyre. 

In table II, the approximate speeds of some of the computed currents are compared 
with data from reference 6. In general, the computed current speeds are somewhat lower 
than observed values. However, a homogeneous model cannot be expected to produce 
quantitatively accurate currents. The low current speeds in the Gulf Stream and the South 
Equatorial Current were attributed to land boundary treatment, grid resolution, and closed 
boundaries, as previously discussed. 

CONCLUDING REMARKS 

Calculations were performed for the near- surface currents of the North Atlantic 
Ocean on a finite-difference grid with 2.5*^ resolution in longitude and latitude and with five 
levels in the vertical direction. The model ocean was homogeneous with a uniform depth 
of 100 m and with a constrained upper surface. The model was driven by published mean 
wind stress data for the Northern Hemisphere summer. The northern (80*^ north) and 
southern (10° south) limits of the model were treated as closed boundaries. 

The equations were integrated for approximately 40 simulation days, at which time 
the momentum equations over most of the ocean were essentially in steady state. The 
computed surface currents agreed for the most part with the observed currents, con- 
sidering the simplicity of the study. Underprediction of current speeds was attributed to 
the lack of density stratification, to the smearing caused by the grid resolution, and to the 
model's treatment of flow along continental boundaries. The development of a strong sub- 
surface countercurrent near the Equator was observed. 

After about 40 simulation days, numerical difficulties at northern latitudes (above 
60° north) led to instabilities in velocity. Minimal applied wind stress in the region of 
difficulty was probably the contributing factor to these numerical difficulties. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
May 10, 1976 
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TABLE I. - MODEL PARAMETERS 


Horizontal grid spacing, deg 2.5 

Number of vertical levels 5 

Depth (uniform), m 100 

Layer thickness (uniform), m 20 

Time step, sec 2600 

Radius of Earth, r, km 6378 

Angular velocity of Earth, co, rad/ sec 7.2919 x 10”^ 

Vertical mixing coefficient, Ky, m^/sec 0.01 

Von Karman constant for subgrid diffusion, Kq 0.4 

Drag coefficient on bottom, Cj) 0.0001 

Temperature (uniform, constant), T, °C 0 

Salinity (uniform, constant), S 35 

Initial velocity (uniform), Uq,, m/sec 0 

Wind stress (constant), Tq, Summer mean 


TABLE n. - COMPUTED AND OBSERVED CURRENT SPEEDS 



Current speed, m/sec 

Current 

Computed, 
present study 

Observed 

reference 

Gulf Stream 

0-5 to 0.75 

0.5 to 2.0 

North Atlantic Current 

0.2 to 0.4 

0.25 to 0.5 

North Equatorial Current 

0.2 to 0.5 

0.25 to 0.4 

Equatorial Countercurrent 

0. 15 to 0.3 

0.25 to 0.5 

South Equatorial Current 

0.2 to 0.55 

0.25 to 1.5 

Eastern North Atlantic Current 

0. 1 to 0. 15 

0.2 to 0.3 
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Figure 4. - Applied wind stress field; averages for Northern Hemisphere summer, interpolated from data 
given in reference 5. Hatched regions indicate artificially closed boundaries. 
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(a) In surface layer (0 to 20 m). 

Figure 5. - Computed currents after approximately 40 simulation days. Hatched regions 

indicate artificially closed boundaries. 









NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON, D.C. 20546 


OFFICIAL BUSINESS 

PENALTY FOR PRIVATE USE $300 SPECIAL FOURTH-CLASS RATE 

BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 



571 OO'l Cl U S 76C6C4 S00903DS 

DEPT OF THE SIB FORCE 

AF WEAPONS laboratory 

ATTN t TECHNICAL LIBRARY (SOL) 

KIRTLAND ARB NH 87117 


TMASTER ; 


If TJndeliverable (Section 158 
Postal Manual) Do Not Return 


^'The aeronautical and space activities of the United States shall be 
conducted so as to contribute to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHJNTICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 

technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS; Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

ADMINISTRATION 


NATIONAL AERONAUTICS AND SPACE 

Washington, D.C. 20546 



